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The increase in the output of papers in the field of mathematical 
biophysics makes it difficult to insure prompt publication without an 
increase in the size of the journal. Therefore, the Bulletin of Mathe- 
matical Biophysics inaugurates the following service: 

Upon acceptance of a paper, the Editor, if necessary, will ask the 
author to shorten the paper to an extent dictated by the requirements 
of a reasonably prompt publication. The shortening should in no 
case reduce the paper to a mere abstract. Such a shortened paper will 
be published within six months or less. 

The unabbreviated original manuscript will be kept on file at 
the editorial office. Any person desiring to avail himself of the com- 
plete manuscript, may obtain promptly a microfilm copy of the latter, 
at the cost of 1¢ per page plus postage, by applying to the Editorial 
Office, 5822 Drexel Avenue, Chicago, Illinois. 

All papers in the Bulletin which have been thus shortened, will 
be marked at the end by the symbol MF, followed by a figure, indicat- 
ing the number of doublespaced typewritten pages of the unabbrevi- 
ated manuscript. 


BULLETIN OF 
MATHEMATICAL BIOPHYSICS 
VOLUME 8, 1941 


A NEURAL MECHANISM FOR DISCRIMINATION IV: 
MONOCULAR DEPTH PERCEPTION 


HENRY E. STANTON 
THE UNIVERSITY OF CHICAGO 


The mechanism for discrimination described by A. S. House- 
holder is applied to the problem of monocular depth perception and 
the Weber ratio is determined. The ratio obtained varies as a small 
positive power of the distance for nearby objects, while for distant 
objects it is indicated that discrimination is relatively impossible. 


In a previous paper, A. S. Householder (1939) has described a 
hypothetical neural mechanism for discrimination, which had the 
property of diverting to different neural pathways responses to stimu- 
li of different intensities. In a later paper (1940) the mechanism was 
applied to the monocular discrimination of lengths which were placed 
at right-angles to the line of sight, and also to binocular depth per- 
ception. (These papers will be referred to hereafter as I and III re- 
spectively.) In both of these cases, the problem was formulated in 
terms of stimuli which, if acting alone and uninhibited, would cause 
a rotation of the eye, and the responses of such stimuli entered the 
mechanism of discrimination. 

The object of this paper is to apply this same mechanism to 
another type of discrimination, that of monocular depth perception. 
In this problem, the change in shape of the crystalline lens during the 
process of accommodation of the eye and the resulting changes in 
muscle tension in the ciliary body are assumed to bear the same re- 
lation to the action of the discrimination mechanism as the rotations 
of the bulbus as discussed in ITI. 

It is, therefore, necessary to make some sort of mathematical 
estimate of the changes in shape of the lens as different points on the 
axis of the eye are fixated in order to determine the changes in muscle 
tension and find the associated Weber function. Only one eye will be 
considered. The problem seems to be naturally divided into two parts, 
(1) the determination of the change in equatorial diameter as a func- 
tion of the radius of curvature of the anterior surface of the lens, and 
(2) the determination of the functional dependence of the radius of 
curvature on the distance of a given point from the eye. 

1. The determination of the change in equatorial diameter as a 
function of the anterior radius of curvature. In Figure 1, an approxi- 
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FIGURE 1 


mate cross section of the lens is shown in the position of partial ac- 
commodation, and this diagram essentially contains most of the geo- 
metrical assumptions made for mathematical simplicity. According 
to W. S. Duke-Elder (1932), in the relaxed position both the anterior 
and posterior surfaces of the lens are spherical. When in the accommo- 
dated position, the posterior surface changes only slightly in curva- 
ture, and will be considered as not changing in this analysis. The 
anterior surface changes very markedly in curvature and moves for- 
ward. The exact nature of the surface which connects the anterior 
spherical cap with the equatorial region is not too definitely specified, 
but indications are that it is practically a cone and has been assumed 
so here. In addition it has been assumed that throughout the range of 
accommodation, points (1) and (2) remain stationary while point 
(83) moves in such a way that the slope of the line (2)-(3), which is 
an element of the cone, remains fixed, the whole process taking place 
at constant volume of the lens. 

In discussing these approximations, they are introduced neces- 
sarily to make the problem sufficiently simple so that it can be solved 
explicitly, and some of them are highly doubtful, although the final 
answer may not be nearly as far from the truth as the original 
assumptions are. In particular, according to most students, point (2) 
does not remain even approximately stationary. This point will be 
common to all possible configurations of the surfaces, whereas it is 
generally assumed that there is no such point, that is, the point cor- 
responding to it would move with the changes. Also the equatorial 
region is probably not circular in cross section, but is frequently 
drawn parabolically. The assumptions regarding the conical surface 
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are probably a greatly over-simplified picture, but the effects on the 
computations will not be greatly altered by more complicating as- 
sumptions in regard to this surface. The important portion of the 
whole problem is the equatorial belt, since it is a result of the deforma- 
tion of this portion which is responsible for the discrimination, and 
this fact must be remembered in any approximation. Unfortunately, 
there does not seem to be very much available literature on the precise 
geometry of the lens in these parts. 

The problem, then, is one of finding the equatorial radius a + r 
as a function of FR, , the radius of curvature of the anterior surface. 
Since it is assumed that the posterior surface does not change in 
shape, the volume V, can be neglected, and the evaluation of the vol- 
ume of the rest of the lens can be considered. For this purpose, the 
lens is assumed to be divided into three parts V,, V2, Vs as shown. 
If K is the slope of the line connecting points (2) and (3), then 


Deere p= KK 
Ue — U3 Le — & 


(1) 


If points (1) and (2) are symmetrical about the y-axis, then from 
the equation of a circle for the curve bounding V;, 


0)? ae 9 (2) 
from which it follows that 
(Ys — @)? + ay? = (Yo — @)? + 29? = 7", (3) 


Another relation is required in order to have a determinate problem 
to connect the coordinates of the points (1) and (2), and various 
assumptions can be made at this point. A very simple one is the as- 
sumption of symmetry described above whereby 


gQ—=— =r 
(4) 
Y¥2=Y=—4, 
ry’ and a’ being the values associated with the relaxed position of the 
lens for the quantities a and 7. 

The coordinates of the point (3) can be determined as functions 
of c and R, quite simply, since the circle is assumed tangent to the 
line segment at this point, 

Y, — R, cos 6, 
(5) 


“x, —-c—F, sin 6, 


where from (1) 
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1 

cos 6 = ——_——_, 

Y aa 

(6) 
: K 
sin 6 = — ——__—_.. 
Vee 
Another relation is obtained by substituting (5) into (1), 

Raga eee (7) 

V1+ EK 


which gives R, as a function of c and the parameter K. 

It is now possible to calculate the volume of the lens having lo- 
cated the important points which are the limits of integration. The 
lens is a solid of revolution and consequently a convenient element of 
volume can be defined as 


dV = 2ay[fi(y) — fel(y) lady, (8) 


where f,(y) and f.(y) are the appropriate equations of the bounding 
curves in Figure 1. 


According to the assumption of constant volume it follows that 
6(Vi+V.+t+ V3) =0, (9) 


and because of the manner of selecting these volumes, V; and V, can 
be considered functions of R, alone and V; a function of a only, and 
so this last equation implies that 


_ a; _ dV, , aVs 


do dhe die eo 
These volumes can be computed from the equation 
V=2aSTA(y) — fe(y)] y dy (11) 


and then differentiated with respect to the proper variable in each 
case and substituted into (10), thus establishing a functional rela- 
tionship between R, and a of the form 


f(R.) + 9(a) =0. (12) 


This expression, theoretically, can be solved for a as a function of 
the radius, but practically this is entirely too cumbersome, and ap- 


proximations will be used. A very simple type of approximation is 
to find the derivative, i.e. 


dR, g(a) (13) 
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and then replace the differentials by the corresponding difference 
values, which gives 


f (Re') 

g(a’) 

The primed values in this expression are assumed to be those 
corresponding to the relaxed position, and this relation will be as- 
sumed to give the dependence of the equatorial diameter of the lens 
upon the radius of curvature of the anterior surface. If the compu- 


tations are carried on in a straightforward manner, the final result 
turns out to be 


dy — A, = — (R2,— Rz,). (14) 


2k?+1 
Zits 1 = SS ee 
2K \/ 1 4 Ke 

which, for simplicity can be written in terms of suitable parameters. 

a=V—-WR.. (16) 


2. The determination of the radius of curvature as a function 
of the object distance. In Figure 2 an approximate cross section of 


FIGURE 2 


the complete optical system of the eye is indicated as far as rays 
which are close to the optical axis are concerned, together with the 
assumed coordinate system. The object of the complete system is lo- 
cated at x, and the image is formed on the retina, the latter distance 
being fixed. The procedure adopted here was to compute the images 
and objects for the successive surfaces using the image of the surface 
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immediately preceding as the object of the surface under considera- 
tion and subsequently finding its image according to the well known 
formula 


Ne ON ees 


: (17) 


Vv wW lie 


The radius of any surface is considered negative if it is convex 
toward the right. The image and object distances are the coordinates 
in the assumed frame and can be either positive or negative. The 
quantities to be substituted for the three surfaces are, 


surface 1 2 3 
object distance, w Oe 2 Or Te 
image distance, v on a Oleg Dae 
where 
P= On ads (18) 


represents the position of the posterior pole of the lens which is as- 
sumed to remain constant or fixed. In these relations, there are only 
three variables, %) which is the independent variable, R, and d,. 

A functional relationship can be established between R. and di 
from the first part of this paper. These quantities when substituted 
into (17) and the functional relation mentioned, which is linear to a 
first approximation, give the radius R, as a function of the single 
variable x, of the form, 

R,=-M+H, (19) 
Xo 
and this may be used with equation (16) to give the value of the 
equatorial diameter a as a function of x) with the form, 


Q 

Cit me (20) 
3. Application of the discrimination mechanism. The paramet- 
ers P and Q can be considered as fixed constants which are deter- 
mined by the geometry of the eye in the relaxed position and they are 
both positive. Since a is also to be considered positive, the formula 
applies only for such values of x, for which the second term remains 
smaller than the first. This point of breakdown of the theory can well 
be anticipated, but this will occur only for points a few millimeters 

in front of the cornea. 


In applying the mechanism described in III, an identical mech- 
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anism can be considered as operative in connection with the ciliary 
muscles which are supposed to give the accommodation. In this case 
the contraction of the muscle which is proportional to the change in 
radius a of the lens is assumed proportional to the activity of center 
C, and this will lead to the expression 


U=v(a — @) (21) 


where u denotes the stimulus intensity, v is a coefficient, and 
where @, would be the value of the radius of the lens in the relaxed 
position, i.e. in the absence of any stimulus. This expression is quite 
analogous to equation (14) in that paper and derived in the same 
manner. Then, if 6 is the Weber ratio for u and A corresponding ra- 
tio for x, , it follows that 


Q 
— (A — Red ed Ses 22 
ut du=v(a, nas atariay>) (22) 
from (72). Hence 
Q Q 
So be se ee ee 23 
Ke ee: | 22) 
But from (11) of III, 
6= pw, (24) 
where a has some value between 1/2 and 1/3. Hence 
5 w= Bue v#(ay— P+ 2+, (25) 
0 


This may be equated to (23) and the equation solved for A a, , which 
turns out to be 
ae nh Ome (Qo =: P ae “4 tee 
Ih oe “ Q : (26) 
aL SF ray a a P == oe 


Assuming that a, corresponds to the relaxed position, it becomes equal 
to P and this equation simplifies to 


14+a 
p %, 


eee eas (27) 
HEL T+ pa’ 


where 
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B 
= y7Qgra, (28) 
WAC 
The presence of the minus sign perhaps requires explanation. 
The formula, in fact, predicts the amount by which an object must be 
moved closer to the subject in order that it may be recognized as 
closer, since moving the object closer results in an increase in the in- 
tensity of the stimulus (deviation of the lens from the relaxed posi- 
tion). 
It is to be noted that for 2, sufficiently small, i.e. for objects suf- 
ficiently close, the J. N. D. will have the approximate value 


Ate = pais (29) 


while as 2% increases, A %) approaches 2, in numerical value, while re- 
maining always numerically less. The interpretation of this is that 
only relatively near objects can be distinguished monocularly as to 
depth—which is, indeed, almost self-evident. 

The writer is indebted to Dr. A. S. Householder for assistance in 
the preparation of this paper, especially concerning the biological as- 
pects. 

This work was aided in part by a grant from the Dr. Wallace 
C. and Clara A. Abbott Memorial Fund of the University of Chicago. 
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SOME PRELIMINARY CONSIDERATIONS CONCERNING 
CONCENTRATION OF OXYGEN IN TISSUE 


INGRAM BLOCH 
THE UNIVERSITY OF CHICAGO 


An approximate expression for the average concentration of 
oxygen in vascular tissue is derived on the basis of certain assump- 
tions. These assumptions are then analyzed, and the inadequacies 
of the approximate derivation are pointed out; possibilities as to the 
method and results of more exact derivations are briefly discussed. 


The problem of what factors influence the concentrations in tis- 
sue of capillary-supplied substances, especially oxygen, was first sug- 
gested to the writer by Dr. Melvin H. Knisely of the University of 
Chicago, in connection with the variations in blood-flow rates during 
malaria. In this paper the writer is attempting merely to outline the 
quantitative aspects of capillary production of metabolites, and to 
derive an approximate expression for the average concentration of 
oxygen in a tissue. He realizes that the treatment here given is by no 
means complete, but hopes that it may, in the future, be made more 
exact and extended to concentrations of other substances than oxygen. 

Oxygen is supplied to respiring cells in tissue by the blood; some 
of it exists in the blood in chemical combination with the hemoglobin 
in erythrocytes, while some is in solution in the plasma. While the 
blood is passing through a capillary, the dissolved oxygen diffuses 
through the thin capillary wall into the surrounding tissue, where it 
is consumed. Meanwhile, the oxidized hemoglobin in the blood cells 
is giving up some of its oxygen to help maintain the concentration of 
dissolved oxygen. 

Obviously the rate of diffusion through the wall of a given cap- 
illary will depend on several factors: 


2, , X= the concentrations, in gm. per cc. of blood, of oxygen in 
arteriolar and venular hemoglobin, respectively. 


¢, = the concentration, in gm. per cc., of dissolved O, in the 
arteriole. 


¢, = the concentration, in gm. per cc., of dissolved O, in the 
capillary just before it empties into the venule. We 
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are not letting c, be venular concentration, because 
in the venule blood from several capillaries is mixed, 
and the effect of the capillary in question is no long- 
er distinguishable. 


h =the permeability of the capillary wall to diffusing O.. 
D=the average tissue diffusion coefficient for O2 . 


—q =the rate of O, consumption in gm. per cc. per sec., in the 
tissue. 


v =the velocity of blood flow in the capillary. 
1 =the length of the capillary from arteriole to venule. 


7) = the radius of the capillary. 
d =the average distance to adjacent capillaries. 


For the present, we simplify matters by making certain assump- 
tions and approximations: We assume that the capillaries in this par- 
ticular tissue are parallel (this condition is approximately realized in 
striated muscle). We further assume that a capillary supplies oxy- 
gen only to the tissue in a right circular cylinder, coaxial with the 
capillary of length / and radius d/2, within which cylinder we are 
neglecting the effect of other capillaries than the one in question. The 
tissue we assume to be homogeneous and isotropic. The oxygen con- 
sumption, —q, we assume to be the same at all points in the region 
of supply of the capillary we are discussing. We are letting the z-axis 
be parallel to the capillary, with z = 0 at the arteriole and z = 1 at 
the venule. The concentration, ¢,(z), of dissolved O, in the capillary, 
is assumed to be constant over any cross section, and to vary linearly 
with z from ¢,; to ¢.; similarly, the concentration, x,(z) of oxygen in 
hemoglobin, is assumed to vary linearly with z from zx, at z = 0 to 
2, at g=1. Finally the external radial concentration gradient is as- 
sumed constant over any plane through the region of supply and per- 
pendicular to the z-axis. We neglect any flow of O, parallel to the 
z-axis in the tissue. We denote by c’(z) the concentration of O. in 
the tissue immediately outside the capillary wall, while é” is the con- 
centration at 7 = d/4, z = 1/2, which we shall call the “average” 
concentration in the tissue. 7 here is the radius vector in cylindrical 
coordinates, measured from the axis of the capillary. 

There are three relations, equations (1), (2) and (3) below 
which are immediately suggested. First, the total efflux of oxyuen 
from the capillary per second may be expressed in terms of arteriolar 
and venular oxygen content of the blood: 
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efflux per sec. =z V 1o2(C, — Cz + 2, — X). 


This quantity must be equal to the total oxygen consumption per sec- 
ond in the region of supply, since we are dealing with a steady state: 
AUT, ¢ + 2, — m2) =— agi] (5)— re. (1) 

A second relation arises from the fact that the efflux per second 
through a small element of the capillary wall between z and z + dz 


may be expressed in terms of the permeability of the capillary wall 
and the internal and external concentrations of dissolved O.: 


efflux per sec. through dz = h[c, (2) PGi ( 2) at T5 OZ. 


This may be set equal to the difference in total oxygen content 
of the blood before and after passing through dz: 
des (Zz) dx, (Z) 

dz = dz 


Pete ca(Z) ~e(@)] de= —| [eared 


or 


(2) 


2h [e.(2) — ¢'(z)]1=— “= ie 2). on, 


Since we are assuming that c,(z) and c'(z) vary linearly between 
the arteriole and the venule, 


déo(Z) Cx" Cy 


dz l 
dio(Z) 2 — Wy 
EOS me ey ei 
Also, 
Cs (Z) == (hy = Aer = Ca) z 


l 


Using these values and the result of (1), we obtain an expres- 
sion for c’(z), the concentration of dissolved O, immediately outside 
the capillary: 


C(e) = Ci a (z)- a s : (z)- a | 


Ure Si a (G5 ae £5) SF 2rh 


Now we set up our third relation, which will lead us to an eXx- 
pression for é", the average O, concentration in the external tissue. 
First, we shall call ¢ the average value of c’(z), which will exist at 
See 
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dC We 


2076" 
Since we assume a linear fall of oxygen concentration along any ra- 
dius vector, the radial gradient of concentration in the plane z = 1/2 
will be given by 


C=, + Se ago) Bet 


where R = d/4. Hence the total flow outward through the surface of 
the cylinder of length J and radius R may be approximated by 
2ate (6 — 6.) D 
R 
Obviously this quantity is equal to the total flow outward through 
the capillary wall, minus the amount consumed in the cylinder: 


2 attD (6 = 6) =a 12 YG Cot a ae) ole eee 
or 


=2na1D(é —€). 


= =" 
= x 


APG Cpa pe) G(T) 
C FS + —___ 


21D 2D 


On using values obtained above for ¢, — Co. + 4%, — X% and ¢’, and 
collecting terms, we obtain 


(3) 


l Re 
c= Cy =i 4 (24 ara is) x q(4 fr 0°) [as : 2D 


Qrh vr? 

This expression is seen to be reasonable, since é” is an increas- 
ing function of r,, h, q (the “rate of production” of oxygen in the 
tissue; q << 0),D,v, and x, — #.; it is a decreasing function of J and 
ane 

Although we have what seems a qualitatively reasonable expres- 
sion for ¢”, it remains to analyze our initial assumptions in order to 
investigate the possibility of a more exact and general approach to 
the problem. 

The first assumption that appears to need refining is that of lin- 
ear change in c,(z) along the capillary. It can be seen that, if there 
were no hemoglobin, all the oxygen being in solution, the concentra- 
tion ¢,(z) would fall exponentially, since, roughly, the rate of fall of 
Co(z) would be proportional to c,(z) at any z. 

The effect of the hemoglobin will, to first approximation, be to 
hold c,(z) nearly constant over quite a wide range of values of total 
oxygen #,(z) + ¢(z) ; this is due to the fact that the hemoglobin con- 
tains many times as much oxygen as does the plasma. 
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However, this approximation will be vitiated by at least two com- 
plicating factors. The first of these is the fact that the equilibrium 
constant of the dissociation of oxyhemoglobin is not fixed; the reac- 
tion is pushed in the direction of dissociation by the presence of free 
carbonic acid. The second complication is the time lag in the dissocia- 
tion, which will prevent the existence of an equilibrium at any point 
along the capillary. These two effects will tend to neutralize each 
other, but just what the final result will be is not known at present. 

Another assumption which will have to be altered is that of con- 
stant radial gradient of external concentration. Poisson’s equation 
shows us that this gradient can never be constant in consuming tissue, 
and we intuitively feel that concentration ought to fall asymptotically 
to zero as 1’, the radius vector in cylindrical coordinates, increases. 

In order to arrive at a more exact idea of the nature of the ex- 
ternal gradient, we must first have consumption, —q, as a function 
of concentration, c. N. Rashevsky (1940) gives the following rela- 
tion for oxygen consumption (the notation has been changed) : 


é=ay -—_—— (4) 


where 
¢ = O, concentration outside a cell 


a, 6 are constants 
y = q/qa*, where —q* = limiting consumption for ¢ = oo. 


If we solve Poisson’s or an equivalent equation for c, having 
first substituted for —q its value as given by equation (4), we should 
get an exact expression for oxygen concentration as a point function 
in consuming tissue. 

Finally, it will be well to discuss another of the concepts em- 
ployed in our approximate derivation of é’. This is the concept of a 
“region of supply’. In our approximation we assumed that the cap- 
illary in question supplied all the oxygen consumed in a region as- 
sumed to be cylindrical; we neglected the effect on concentration in 
this region of all other capillaries than the one under investigation. 
Although it is true that, if there are n capillaries supplying a volume 
V of tissue, the average capillary may be assumed to supply a volume 
V/n, and although it is further true that through any closed surface S 
constructed around any such volume V/n containing one capillary the 
net diffusion is zero, it is not true that the capillaries outside S have 
no effect on concentration inside S. It is obvious, for instance, that 
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if one of the external capillaries stops flowing, S will, on the average, 
have to expand to include the volume V/n — 1; hence, concentration 
within S will be reduced. 

The problem of determining the exact effect of distant capilla- 
ries may be approached in either of two ways: by summing, or inte- 
grating, the effect of a single capillary over all the capillaries in the 
tissue, or by finding the total diffusion out through S of products of the 
capillary within S. This quantity will be equal in magnitude to the 
diffusion inward from all other capillaries, the average value of which 
diffusion will thus be determined. 

In concluding, the writer wishes to express his appreciation of 
the very helpful assistance and counsel given him by Dr. A. S. House- 
holder. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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NOTE ON A POSSIBLE APPLICATION OF SOME CONCEPTS OF 
TOPOLOGY TO ASYMMETRIC ORGANIZATION 
OF PROTOPLASM 
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A suggestion is made for the possible use of topology in the 
theory of asymmetric constituents of protoplasm. 


At the present time we are not provided with the exact knowledge 
of the molecular skeleton of living matter, which could be placed at 
the bottom of any elaborate mathematical theory of structure and 
activity of living systems. The only starting point in our possession is 
presented by observations on optical activity of living matter. It 
appears reasonable in this connection to use some concepts of a 
specialized branch of mathematics, namely topology, in order to ex- 
tract from the empirical records concerning optical activity of proto- 
plasm some more general information concerning the properties of 
space in living systems. Topology is the science of general properties 
of space, and some of its generalizations have repeatedly been found 
useful in all those fields of natural knowledge that deal with the 
properties of complex spatial figures. 

The empirical data of biology show that the molecules of various 
substances participating in the fundamental metabolic processes 
always possess the same (left) steric configuration (for a summary 
see Gause, 1941). It is only in leaving these fundamental processes 
that the organic molecules can sometimes change their configuration 
into an inverse, the right one. Such an inversion of configuration is 
often observed in pigments, excretory and some other secondary sub- 
stances. In other words, the primary physiological processes involve 
molecules of the same left orientation. These molecules form an equi- 
oriented manifold, in the topological sense. However, the equiorienta- 
tion does not extend upon all the manifolds of molecules in living 
systems, but is restricted to primary constituents only. In the sec- 
ondary substances, dissociated from the primary complex, the equi- 
orientation does no longer exist. 

The equiorientation of a manifold belongs to one of the basic 
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concepts of topology. All the manifolds of three-dimensional figures 
(molecules included) are divided by topology into two categories: 
first, the orientable ones, which can be oriented throughout in a single 
specific manner (equiorientated); and, second, the non-orientable 
ones, which cannot be oriented in the sense just given. 

In this way from the view point of topology the diversities be- 
tween the primary and secondary constituents of living matter are 
very profound, as far as these belong to two different categories of 
topological structures. The primary constituents can be characterized 
as the orientable manifolds, whereas the secondary constituents repre- 
sent non-orientable manifolds. 

Having now passed from biological to topological language, we 
can analyze a little further the diversities between the orientable and 
non-orientable manifolds of molecules. In the equioriented manifold 
subject to various transformations any inversion of orientation of 
a component figure is excluded. Topology teaches us that such an ex- 
clusion requires some special conditions. It is required that such trans- 
formations should be “forbidden” which by purely geometric reasons 
inevitably lead to an inversion of orientation of a figure. A relatively 
simple geometric example of such a “forbidden” operation of inversion 
is provided by the so-called Mobius band (Seifert and Threlfall, 1934). 
When a circle rotating clockwise is passed through such a band and 
then placed against its previous image, it will now be found rotating 
in a counter-clockwise direction. 

In other words, the equiorientation of the manifold of molecules 
in primary constituents of protoplasm is attained by means of for- 
bidding some specific intermediate topological transformations of 
these molecules. Such a forbidding is not operative for the secondary 
constituents of protoplasm. 
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Electric impedance measurements on systems consisting of non- 
uniform elements (e.g. nerve bundles, cell suspensions, imperfect 
dielectrics) can be interpreted if the impedance characteristics of 
the individuals and their impedance distribution are known. Con- 
versely, given the observed overall impedance, z(p), the impedance 
distribution of the individuals is not uniquely determined in both 
phase angle and static capacitance. If all individual phase angles 
are equal, the distribution, D(r), is the solution of Stieltjes’ integral 
equation 


"co D(r)dr 
Fiew)1= | sion tas 
S 1+ pr 


F being a known function of the geometry of the system. It is shown 
that the phase angle of the equivalent suspension approaches that of 
the original suspension as the dispersion of D(r) decreases. 


The majority of impedance measurements on biological materials, 
have, until quite recently, (Cole and Curtis, 1938b) been made on 
multifibered preparations (e.g., nerve or muscle), or multicellular 
suspensions (e.g., marine eggs or erythrocytes) (Cole and Curtis, 
1936; hereinafter abbreviated C 1936). In general, the aim of such 
measurements has been to determine the impedance characteristics 
of the individual units comprising the system from the overall re- 
sults found for the entire population. In the case of certain marine 
eggs (Hipponoé and Arbacia, e.g.) where the population is nearly 
uniform, Maxwell’s calculations for the impedance of suspensions of 
spheres have given results which are entirely consistent with imped- 
ance measurements on single eggs. On the other hand, nerve and 
muscle, which consist of parallel fibers showing wide variation in 
diameters, membrane capacities, and phase angles, give rise to overall 
“nolarization” impedance characteristics which do not agree, insofar 
as the comparison has been made, with experiments on individual 
axons. The suggestion has been made (C 1936) that the appearance 
of polarization impedances in biological materials may be attribut- 
able entirely to the statistical variation of the electrical parameters 
(other than phase angle) of the individuals of the population. This 
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interpretation is not completely adequate in the only case which has 
been investigated, viz., the fin nerve of the squid; for, as Curtis and 
Cole have pointed out, (1938b) individual giant axons exhibit a po- 
larization impedance with an average phase angle of about 76°, 
whereas the phase angle for an equivalent uniform nerve lies between 
40° and 50°. It thus appears that the observed polarization imped- 
ance in biological systems is partly due to polarization impedance of 
the individual members, and partly due to a statistical distribution of 
sizes, static capacitances, and phase angles in the preparation. 

The statistical distribution of the individuals in a system whose 
overall impedance characteristics are known, is related to the imped- 
ance of the equivalent homogeneous system by an integral equation, 
which, in certain simple cases, reduces to an equation recently dis- 
cussed by R. Fuoss and J. Kirkwood (1941) in connection with the 
distribution of relaxation times in imperfect dielectrics. However, 
since the integral equation is an instance of one treated by Stieltjes 
in connection with his moment problem, it may be of some mathe- 
matical interest to present-an alternative solution as an application 
of the Stieltjes theory (Perron, 1929). 


The Distribution Problem 


The present considerations will be confined to the case mentioned 
in C 1936: the impedance of a statistical distribution of parallel cylin- 
ders, although the final results are applicable to spheres with very 
little change. Let us suppose that the system consists of three phase 
(external, membrane, internal) parallel cylinders each of whose mem- 
branes has specific impedance z(7w)-”, where w/(2x) is the fre- 
quency of the measuring E.M.F., m 2/2 is the membrane phase angle, 
and i= ./—1. Then by a simple extension of some considerations in 
C 1936, the specific impedance, z, of the group of fibers whose vol- 
ume concentration, G(z, m), lies between G(z, m) and G(z + dz, 
m + dm), is determined by 


Det PAN in aes Hopi @) ane 
a( : \=G@,m) pie citi NIE Ss LG 
@ ate Ty fags +e Hey + z(4 ay ; (1) 


where r, and 7, are the specific internal and external] resistances, re- 
spectively, and G is normalized so that 


[°-[ G@,m) amdz = 6, (2) 


G, being the volume concentration of the equivalent uniform suspen- 
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sion. The specific impedance, Z , of an equivalent homogeneous sus- 
pension is given by 

foe eT Tet z; 

fae ee Re ae eer aoe Q 

Z+r, a +7, + 2s” (3) 
where z; is the equivalent membrane specific impedance. If we equate 


the right member of (8) and the integrated right side of (1), we find, 
after subtracting G, from each side, 


=| ee G(z,m) dmdz 

1 Pear G.[1+ za) 77) es 
whereW=7n4+7. 

¢ The function 2;(w) is known empirically to have the form 
23 (i w)-*, z, and f being constants. However, it is easy to show that, in 
general, (4) does not determine G(z, m) uniquely. For, suppose the 
distribution in z and the distribution in m are uncorrelated. Then 
(Kenney, 1939, p. 64) G(z, m) must be expressible as the product 


of two functions, h(z)g(m). Substituting this product for G(z, ™) 
in (4), and integrating with respect to m, we find 


1 7 *00 — tes == 
moe K(o,2)h(2) dz, (5) 


where the kernel, K (@ oe is defined by 


g(m)dm 


J ate oe 
o[1 + (io) -"/W] 


Kioes) = 
Since the solution, h(z), of the integral equation (5) depends on 
K(@ ,z), which in turn depends on g(m), we can choose g(m) arbi- 
trarily and get a different h(z) for each admissible choice. This means 
that data on z;(w) alone will not determine G(zZ, m) ; it is also neces- 
sary to know the distribution in one of the variables, z or m, in order 
to find the distribution in the other. 

In the absence of any data on g(m), we shall make the simplest 
assumption (which is certainly incorrect for nerve, but is probably 
correct for marine eggs), viz., g(m) = 6(m— mp) where 6 is the 
Dirac delta function. In this case equation (5) may be reduced to a 
canonical form 


ioe (i D(z) dz (7) 
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a= f/m, 
p= (ta) (23/W) 2, 
r= Wiv/e-0z/2,1/2, 
D(z) = (1/Go) 23/4 W/™ h(z). 


The physical significance of the distribution variable z depends on 
the impedance characteristics of the individual membranes. Thus, for 


(3) 


elements whose membranes are pure static capacitances, z = (Cua)" 
and m= 1, where Cy is the membrane capacitance per unit area and 
a is the fiber radius. If, furthermore, the specific properties of the 
fibers are all alike so that Cy is the same for all fibers, the distribution 
will be a simple size distribution. However, the assumption m, = 1 
and z == (Cya)- is somewhat irrelevant for nerve fibers in the light 
of Cole and Curtis’ single fiber measurements; it seems to be justified 
only in the case of marine egg suspensions (Cole and Curtis, 1938a). 


Solution 
The substitutions 


$(1) = [ Doar, pace (9) 
reduce (4) to 


cet .0 (7) 

ee ages SaLAGES (10) 

Ih ae (ce C+ 

which is the form associated with Stieltjes’ moment problem. We em- 

ploy Stieltjes’ method of solution as given by Perron (1929, p. 372). 
This solution ¢(&), where continuous, is given by 


cae Oem Wake Set! a 

(€) = 29(0) =lim # — | —— 
o  ( ) ae i cae 1 + ce ? (11) 
where the value of the function FR is the real part of its argument; 
from (11), this is —1/z times the coefficient of 7 in the definite inte- 


gral. For any 7 # 0 the integral has the value 

+in 1 ) a 
Se Lee Ea 1) : 

eae Tae (2h 

and the limit of the imaginary part of this logarithm as 1 approaches 


zero from the positive side is needed for (7). Since the argument of 
the logarithm can be written as 


Riba (1 + ¢*) 
a 


(12) 


\ 


Pet (01-62) 
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where P is real and #, and ¢, are the amplitudes of the numerator and 
denominator, respectively, of the fraction, our limit becomes, upon 
taking logarithms, 


: 1 
lim —— (6, — 62). 
Qa 


n=+0 I 
The values of 6, and 6, are, from (9), 
yn“ sin a 2/2 


G7 —- arc tan ——___—__ 
; 1 + 9*-cosia.2/2- 


p? sin a 6 yey 
6, = are tan —_—__—_—_,, 
Lp? 208,.0)9 
where 
p—& + 7?, tan 6=-—7/é (14) 


and @ lies in the second quadrant. The quantities 7% and p* are positive 
and real. 

As 7 approaches zero, tan 6, approaches zero, and hence 6, ap- 
proaches zero or some integral multiple of 7. We may, however, ig- 
nore additive constants [absorbing these into ¢(0)], and regard the 
limit of 6, as zero. The angle @ has the limit a, and p? the limit &, 
whence 


: &*sinaa 
lim 6, = are tan ———__—__—_—, 
1=+0 1 + ¢* cos an 


and therefore, where ¢ is continuous [cf. equation (11)], 


&*sinaz 


1+ &*cosan’ oe 


i 
$(é) =¢(0) + =; are tan 
In case a = 1, the arc tangent has the value zero or an integral 
multiple of 2; this implies that ¢(&) is constant except at points of 


discontinuity. The original equation (10) becomes 


(16) 


1 eae (53 d p (7) 

Leg } (Pee 
and if the discontinuities of ¢ occur at the points 7;, then the integral 
has the value 


Rep Aal vain traits 0) dtr 0) 
C+ C+ ; 


0 


This will equal the left side of (16) if and only if there is a single 
discontinuity at 7 = 1, with 
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6(L +10) = 6(1 = 0) = (17) 
If ¢(0) =0, then (cf. Fig. 1) 

6(é) =0, 0 aly 

o(é)=1, bia Ac Ge 


Thus D(r), which is the derivative of ¢(7), is zero everywhere except 
at +r = 1, where it becomes infinite. The function D(r) in this case is 
identical with Dirac’s 6-function. 

The physical significance of the discontinuous integrated distribu- 
tion for a = 1 is the following. If the phase angle of the individual 
elements is the same as the phase angle of the equivalent suspension 
(6 =m,), so that a= 1, then, by the previous paragraph, the standard 
deviation of the distribution is zero. This is fairly well verified in the 
case of the marine eggs. Cole and Curtis (1938a) find phase angles in 
single marine eggs corresponding to values of m very near 1, while 
suspensions of marine eggs in which “the standard deviation of cell 
diameter is of the order of one per cent... give a value of 6 = 1.0 
within experimental error” (C 1936). 

For 0 < a < 1, we obtain D(r) by differentiating (15) (the as- 
sumption that ¢(¢) is continuous being justified on physical grounds) : 

1 7 sin an ; 
DG) 5 a oe COR a ua) 
and this is the required distribution. The validity of this formula in 
the special case a = 1/2 is readily checked by direct integration of 
(7). 

The distribution (19) is most conveniently visualized by making 
the substitution s = log 7. Then the logarithmic distribution function, 
F'(s), defined by 


F'(s)ds = D[r(s)] dz(s) 
is, from (19), 


1 sin az 
HAS} 2a coshas +cosan’ 
which is exactly the distribution obtained by Cole and Cole (1941). 
Graphs of ¢(7), D(r), and F'(s) for several values of « between 0 and 
1, are given in the accompanying figure. It will be noticed that as a 
approaches unity, the dispersion of the distribution decreases; 1.e., as 
the uniformity of the suspension increases, the phase angle of its 
members approaches that of the equivalent homogeneous suspension. 
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This result is intuitive, since the more uniform a suspension, the more 
it resembles its equivalent homogeneous suspension. 

The writers wish to thank Professor H. S. Wall of Northwestern 
University for having pointed out to them Perron’s treatment of the 
Stieltjes problem. It is to a remark made by Professor Wall in a lec- 
ture on continued fractions (1941), and to a comment by Professor 
K. S. Cole that the origin of this paper is due. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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Continuing the investigation previously introduced, it is shown 
here that when the product of the activity parameters of the cir- 
cuit is not exceeded by unity (algebraically) a steady state is not 
possible in which all fibers of the circuit are active, whereas when 
this product is exceeded by unity, any stimulus pattern which is con- 
sistent with such a state of complete activity is inconsistent with 
any state of partial activity of the circuit. 


In two previous papers (Householder, 1941; these will be de- 
noted by I and II) the general problem of the nature of the steady- 
state activity in a network of neurons was introduced. The case of a 
constant stimulus only was considered, and it was assumed that a 
single fiber, when acted upon by a constant stimulus which exceeds a 
certain threshold, will then react with an intensity which is asymp- 
totically proportional to the amount by which the constant stimulus 
exceeds this threshold. The threshold and the constant of proportion- 
ality, which was called the “activity parameter,” are characteristic 
of a given fiber. Under these assumptions, the case of a simple circuit 
of n fibers was discussed, and it was shown that to prescribe the 
intensities of the stimuli which are applied from without to the sev- 
eral fibers does not suffice in all cases to determine which of the 
fibers of the circuit will respond, in case there are inhibitory fibers 
present. The results were, however, incomplete, inasmuch as it was 
assumed that there was always at least one fiber of the circuit which 
was inactive. It is the purpose of this note to supplement these results 
by removing this restriction. 

The notation and terminology is the same as that employed in 
I and II, and will not be defined again here. Let 


Aart = Ode <7 dy, , (1) 


where the a; are the activity parameters of the fibers of the circuit. 
Subject to the restrictions mentioned, it was found previously that 
when A> 0, and also when the circuit consists of excitatory fibers 
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only (every a; > 0), then the external stimuli (“stimulus pattern’) 
do determine uniquely the activity of the entire circuit and its amount 
for each fiber (the activity pattern and the excitation pattern, respec- 
tively), whereas in the contrary case, there always exists a stimulus 
pattern with which each of two distinct activity patterns is consistent. 
Here we shall show first that when A < 0 it is not possible for all the 
fibers to be acting at once in a steady state, and second that when A 
> 0, no stimulus pattern which is consistent with a state of complete 
activity can be also consistent with a steady state in which only some, 
but not all of the fibers are acting. 

The first result in the case all a; > 0 has already been obtained 
in I, and simply means that for large stimuli and for such circuits, 
the activity parameter must be regarded as variable, the assumption 
of constancy being an inadequate approximation to the true situation. 
Where inhibitory fibers are present, however, the equations indicate 
that as the stimuli are increased, whichever inhibitory fiber becomes 
active first, this one will tend to dominate the circuit and render 
succeeding fibers inactive, so that in no case (when A < Q) is it pos- 
sible for all fibers to act simultaneously. To see this, consider any 
SPo which is supposed to be such that a steady state of complete ac- 
tivity can result. Then there must exist some mode of gradual appli- 
cation of these stimuli, perhaps a little at a time to one and then an- 
other of the synapses, until all fibers are active. In fact, it should be 
possible to do this slowly, no increment being made to the stimulus at 
any synapse until the steady state has been reached as a result of the 
preceding ones. Then at some stage, all fibers but one can be made 
active, and the final increment will be applied to this fiber to render 
it active, after which all fibers will be acting. But, the final step will 
therefore be to make o;‘”) positive where it was previously negative or 
zero (zt being the number of the fiber). But if this is done, and if all 
fibers thereby become active, then in the steady state the resulting 
values of the y; at the several synapses are the solutions of equations 
(7) of I, where every a; = a; , and this solution is 


Yio 7A, (2) 


as given in (13) of I. But o; > 0, A < 0 by hypothesis, wherefore 
yi < 0 (or else y; = 0), whence the fiber 7 cannot remain in activity 
if all the other fibers remain active. Thus, putting fiber (i) into ac- 
tivity necessarily stops at least one of the other, hitherto acting fibers. 

Now suppose A > 0 so that it is possible to have a steady state 
of complete activity, and let « be an SP with which this state of com- 
plete activity is consistent. The AP is defined by a; = a; and the EP 
is defined by the y; given in (2). Let this same SPo be consistent with 
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some AP £ # a, i.e. an AP which has at least one zero. Let the cor- 


responding HP be z;. Let the zeroes of £ be located at i, laa te 
where 
EES SSeS ce a eS a (3) 
Let 
te Ge —— IP ty, til Oia == toy tra + Wp=%y. (4) 
Also let 
Ai, jaa = Wi Wins ++ A; when t<j<n, 
Ais, jay = Gj Wig ++ Ay +++ 4; When j<i<n. O) 
Then if h does not coincide with any of the7z,, we have 
Yar = Onsr + Yn Zhv1 = One + Anzn, 
whence 
Ya — Zr =A (Ya — Zn) , 
and by induction 
= Ying — Sng = An, mi (Yn — &n) , (6) 


provided the sequence h,h+1,---,4+7-—1 does not include any otf 
the z7,. On the other hand, if 6,—0, while the sequence h + 1,---, 
h + 7 — 1 does not include any other 7,, then z, must be omitted 
from equation (6) and we have 


Ynrji Sng =An pj Yn - (7) 


These equations are not restricted to the case when every a; # 0, 
but are more general. In fact we have the 

LEMMA. If o; is an SP consistent with an AP a; if 7; is an SP 
consistent with an AP §; if 


Onis Onto *** Onej—1 Pave Bresso Bide Pipes #0 > 


and if 
Ont. — Thti 9 °** » Oh+j-1 — Th+j-1 » 
then equation (6) holds provided also mf, #0, but (7) holds when 
Bi, = 0. 
If, now, we seth =i,, and j —o@,., then we have from (7) that 
Y eaees an aia (8) 


Since Chie UZ 


tre 


£0; 4, > 90, wWfollows from (8) that 
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A B(s: (9) 


The equations (8) also imply the inequalities 
Ue en ae (10) 


and from this system of » inequalities we easily obtain 


Y,. s Agus Y,. 9 
and therefore 
A<0, 


contrary to hypothesis. 

It follows, therefore, that when A > 0, every state of partial or 
complete activity of the circuit is possible and is uniquely determined 
by the applied stimulus pattern, whereas for 4 < 0 a steady state of 
complete activity is never possible. Consequently the restriction which 
was repeatedly imposed in II, that each AP considered shall have at 
least one zero, was either unnecessary (when A > 0) or automatically 
fulfilled (when A < 0). 

This work was aided in part by a grant from the Dr Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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On the basis of Rashevsky’s nerve excitation equations, an ex- 
pression is derived for the distribution of response times attributing 
the variation to the fluctuations in threshold. The resulting equation 
is compared with available data and agreement is found. 


The suggestion has been made that the variation in response 
times of a nerve fiber is due to the fluctuations of threshold (E. A. 
Blair and J. Erlanger, 1933, 1936; C. Pecher, 1939). A discussion of 
possible sources of the threshold variations is given by C. Pecher 
(19389). Of particular interest is the suggestion that the principal 
source of the fluctuations is the variation in the fiber as would be ex- 
pected from the kinetic theory. Calculations showed this to be very 
plausible. Various experiments showed that the fluctuations were 
random, and the fluctuations in two adjacent fibers were independent. 

It was shown by H. Jasper and T. Perkins (1932), that for nerve 
muscle preparations the cumulative response curve plotted against 
time could be closely approximated by an ogive. This would imply 
that the derivative curve, or the curve of the distribution of response 
with time, would be a normal curve. However, this is not quite so. It 
is the purpose here to derive this relation more exactly and on the ba- 
sis of equations which have considerable experimental support. 

Let us assume that the variation in response times is due entirely 
to the variation in threshold and using N. Rashevsky’s form of the 
two-factor excitation theory (Rashevsky, 1938, chap. xvii), derive 
the expression for the distribution of response times. If the varia- 
tions in threshold were due to fluctuations in ion concentrations as 
suggested, the latter would be randomly distributed, and since the 
numbers are fairly large, one obtains for the distribution of the thres- 
hold, the well known error curve. It is assumed in the derivation that 
the numbers dealt with are large, which they are, and that the varia- 
tions are small compared with the average value. In terms of the 
threshold distribution this means that the standard error must be 
small compared to the mean threshold. But, this is true for the nerve 
as C. Pecher (1939) obtains, for the ratio, values of about one or two 
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per cent. Thus, we may assume a normal error curve. 
For constant current excitation, N. Rashevsky’s equations are 
(Rashevsky, 1988, p. 164) 


ee (1 —.e""), (1) 


cla ee MI 
J = Jo = (1 ae Gwe) ? (2) 
m 


where the condition for excitation is e =>j. The time ¢ is measured 
from the application of the stimulating current I. 


Let 
eran | 
bm SP 
so that there is normal accommodation, (Hill, 1935). Then 
; : KI 
9 == —" Ga — E5) ray (em teri) a (4) 


The value of ge and 7,, or actually their difference, is considered to be 
constant. But let us now consider that at any particular instant it 
may have some value & — jo — & So that e, — Jo is the average value 
and € is a fluctuation from the average. Then, instead of (4) we have 


: , I 
Banh fection ©) Youn.) ete Re 4 oer 8 (5) 


If p(é) d é is the probability that € has a value between &€ and é + 
d&, then as é is normally distributed, 


1 a 
p(¢) =——— e*™*_, 

V 21 o e 
where o is the standard deviation and equals to 1.67 times the prob- 
able error. The condition for excitation to take place is e > 7, or be- 
cause of (5) we may also write for the condition that excitation occur, 


KI 
eee aie LG a BOT) a ete lie (7) 


If we define P(I,t) as the probability that a response will occur for 
a fixed time ¢ when the stimulus intensity is J, then by the definition 
of P and from equation (7), 
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& 
PUI, t) = | p(é) dé, (8) 
where Ds 
KI : 
(1, t) =—— (e™' — e) — (je). (9) 


If ¢ is constant, then &, is linear with J. But because ¢ is a func- 
tion of 7, €, in (9) is not linear with J. But if one recalls that for 
stimuli near the threshold (it is only these in which we are interested 
for the present) a response will not be apt to occur until &, is very 
nearly zero, that is, until it is a few per cent of the initial value 
(Jo — &0).-From the fact that the stimulus is near threshold, the factor 
emt — e-*t of the expression (9) is near its maximum value and thus 
it is approximately a constant. Hence for values of J near the thres- 
hold (defined as that value of J for which the proportion of responses 
is one-half), &, is essentially independent of t. The unit in which the 
stimulus is expressed is immaterial as long as the threshold and stand- 
ard deviation are measured in the same unit. Thus, defining 


KI 


Aes re (e% ae. err) (10) 
Po =Jo — &o (11) 
then equation (8) becomes 
»J!-I'q 
PL) = paé. (12) 


-o 


A comparison of this expression with data by C. Pecher on a 
single fiber from a frog sciatic, is shown in Figure 1, where the ex- 
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perimental values are solid circles, and the curve is that of equation 
(12) with J’, = 100, and o = 1.667, or 
6 
V 2a 
The agreement further emphasizes the validity of the assumption of 
a normal distribution of € in equation (6). 

We now proceed to the derivation of the expression for the dis- 
tribution of response times for a constant J. If a fiber being stimu- 
lated by: an intensity J , has not responded at a time ¢, the probabil- 
ity that a response will occur between ¢ and ¢ + dt is given by Pdt 
where P is given by equation (8), and where &,(¢) is given by equa- 
tion (9). In order that a response should take place between ¢ and 
t+ dt, it is necessary that no response shall have occurred between 
zero and ¢ (or for that matter for a sufficiently long time before t= 0, 
but no response occurs since J is then zero). Let the time from zero 
to t be divided into n equal intervals, 4d t. If q(J, t) At is the 
probabilty that a response occurs during the interval ¢ to t+ At, 
when the stimulus intensity is], gq A t is given by the product of the 
probabilities that no responses occur in the preceding intervals and 
the probability that response does occur in the interval t tot + At, or 


& 


» (1’/-100) -.86—~— 
| @ 20es (13) 


-O 


(1) 


qAt=(1—P,At)(1—P, At)---(1— Pea 4t)P, At, (14) 


where P; is the probability P(UJ, t%) = P(1,iA t) that a response 
will take place between ¢; and ¢; + A ¢t = ti, if no response has oc- 
curred before ¢;. Equation (14) may be written 


log g=> log(1 — P; At) + log P (15) 
i=0 


where q and P are taken att. Since A ¢t may be taken as small as de- 
sired, we may expand each logarithmic term and use only the first 
term of the expansion in each case, or 


logq=—-DPiAt+logP, 


i=0 
‘ (16) 
—— { Pdz+logP, 
0 
so that 


ie P(1,tT)dt 
—Jo 


qd tyes P Ud, te , (17) 
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where P(I, t) is given by equation (8) and &,(t) of the latter equa- 
tion is given by equation (9) for any intensity of the stimulus J. From 
equation (17) we see that the mean response time must decrease as 
the intensity I increases, a fact which is well known. The spread in 
the response times must get smaller as the intensity is increased. This 
is shown in the experiments by C. Pecher (1939). Whether or not 
the agreement here is also quantitative remains to be seen. 

We next inquire as to the relation between g and t. C. Pecher 
(1939) gives two such curves in which the intensity of stimulation is 
near threshold. As the exact value is not explicitly stated, it was as- 
sumed that the value in each case was exactly threshold and the area 
under the curves in each case reduced to one half. This led to the 
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points shown plotted in Figure 2. Since the values of the probable er- 
rors relative to the threshold were not given, these were left to be 
found by trial and error in fitting the respective curves to equation 
(17). Similarly the values of k and m were not given, though the 
value of k can be obtained to a first approximation without difficulty. 
The ratio k/m has some value from 10-100 (Katz, 1939). 

For the left curve, the following values were taken as parame- 
ters: the ratio of the standard deviation, «, to the threshold (given 
by I in equation (9) for €, = 0 and such a value of ¢ that em’ — e™ 
is a maximum) was 2.16%, or 1.3% probable error; k — .95 msec., 
and m = .0633 msec. or k/m = 15; thus 


— 130 [.7692 + e--95t as Fan | ‘ 
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For the right curve the following values of the parameters were used: 
the ratio of the standard deviation to threshold was .88%, or in terms 
of the probable error .538%: k = 1.35 msecs., m = .027 msec. or 
k/m = 50; or 


E, = 105.4[.9048 + e188 — e-oert] 


The curves are not particularly sensitive with regard to m. A 
fairly large change in m, as from 1/5 k to 1/15 & showed no appreci- 
able effect on the curve, except that in the former case o had to be 
taken as 3.1% whereas in the latter o = 2.16% in terms of the thres- 
hold. The former is a rather high value, the latter is about average 
(C. Pecher, 1939). Thus with the value of c known, there remains 
not much variation possible for m. But as k and m can be indepen- 
dently determined, it is clear that in principle all the parameters in 
equation (17) can be determined independently, and with these given, 
one can plot the expected distribution. Whether or not the actual val- 
ues agree is then readily tested. Under these conditions, agreement 
would strongly suggest that the variations in response times is due 
almost entirely to the variations in the threshold. With a knowledge 
of the latter and sufficient additional information, the distribution of 
response times is predicted quantitatively. 

It is of interest to note that regardless of the value assumed for 
the per cent standard deviation no fit can be obtained if m is taken too 
small. In the first curve, a value of m equal to k/50 is already too 
small, in the second, k/100 is probably too small. Thus to obtain any 
agreement at all, it is essential that a two factor theory rather than 
a one factor theory be assumed. It has been suggested that the second 
factor does not occur under physiological conditions or if the fiber is 
amply perfused (Parrack, 1940). If so, we should expect the distri- 
bution of the response times under the latter conditions to show a 
slower fall for larger values of the response time. The experimental 
curve should be fit by equations (17), (8) and (9) where m is set 
equal to zero. In this case, it should be noted that the threshold as de- 
fined by jo — & is not equivalent to the definitions in terms of a fifty 
per aa response, and one cannot assume them to be approximately 
equal. 

Perhaps it should be.emphasized that the validity of equation 
(17) does not depend upon the initial equations in ¢ and 7. If one is 
given the distribution of thresholds of a particular nerve under a set 
of conditions, and also the empirical excitatory state as a function of 
time, then the operations indicated by equation (17) lead one to the 
distribution of response times. No unknown parameters are involved. 
If the conditions were maintained constant throughout the experi- 
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ments, agreement should be very significant, especially if various 
types of fibers were used and if the initial conditions were varied 
widely. 

An expression for the distribution of response times for a single 
fiber has been obtained under the assumption that this variation is 
due only to the known variation in excitability. The predictions have 
been compared with available experimental data in which not all the 
required data were available. If the data were sufficient the resulting 
expression would contain no arbitrary parameter. In the cases tested, 
the values assumed for the parameters not given are consistent with 
other data. To the extent of this agreement, it would seem that the 
variation in the time required to set up an impulse is largely due to 
the known variations in threshold. 
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APPLICATION OF THE GENERAL FLUID CIRCUIT THEORY 
TO THE CHLORIDE-BROMIDE-DEUTERIUM 
OXIDE EXPERIMENT 


H. C. PETERS 
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The general fluid circuit theory of active chloride absorption is 
applied to an experiment in which chloride, bromide, and deuterium 
oxide were simultaneously absorbed from the lower ileum, chloride 
only moving against a concentration gradient. The conclusion is 
reached that the constant for active absorption of deuterium oxide 
is equal, within reasonable limits of error, to the constant for active 
chloride absorption in accordance with the assumption that the ab- 
sorption of water carries chloride out of the intestinal lumen without 
changing its concentration. 


The original simple form of the fluid circuit theory (Ingraham, 
Peters, and Visscher, 1938) assumes that the absorption of water at 
the rate FR, carries chloride out of the intestinal contents without 
changing its concentration from that in the lumen while a fluid free 
of chloride moves from the blood into the intestinal lumen. If this 
theory is correct and deuterium oxide is a perfect indicator of active 
water movement, then R, calculated from deuterium oxide absorption 
(Rov,o) Should equal R, calculated from chloride absorption (Roc) for 
an experiment in which deuterium oxide and chloride are simultane- 
ously absorbed, chloride moving against a concentration gradient. Pe- 
ters and Visscher (1939) found that R.p,o is greater than Ry: and at- 
tributed this discrepancy to diffusion of deuterium oxide through pores 
other than those involved in active water movement and to the pres- 
ence of chloride in the fluid passing into the intestinal lumen. In an 
attempt to evaluate the importance of diffusion in deuterium oxide 
absorption they carried out an experiment in which chloride, bromide, 
and deuterium oxide were absorbed, chloride only moving against a 
concentration gradient. The general fluid circuit theory (Peters, 
1940) which includes diffusion and secretion of chloride and osmosis, 
suggests an analysis of this experiment which furnishes new evidence 
for the theory. 

It is assumed that bromide and deuterium oxide are each ab- 
sorbed by two first order processes proceeding side by side so that 
for either substance 
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dI/dt =C(dV/dt) + V(dC/dt) =—(atd)C, (1) 
where I igs the amount of the substance in the intestinal lumen, C its 
concentration, V volume, and ¢ time. It is assumed that a is the con- 
stant for an active process and d is the constant for transfer by dif- 
fusion. Assuming further that 


Ve AV gel JE (2) 
where D is the rate of volume decrease, we obtain by integration from 
equations (1) 

CV /CoVo— (V/V e) 2; (3) 
where V is given by equation (2) and C, and V, are initial values for 
C and V. 

From the original fluid circuit theory Peters and Visscher have 
derived the equation 
GV/G5V 5 = (V/Ve) F?; (4) 


where V is given by equation (2). Since equations (3) and (4) have 
the same form, R,, calculated for bromide or deuterium oxide, has 
the significance indicated by the equation 


R,=a+d. (5) 
For bromide and deuterium oxide respectively we have 
Fopr = Apr + Apr , (6) 
Reno = as ae ; (7) 
Using equations (4) and (2) Peters and Visscher (1939) found 
Hone = 3.9hCC./min., (8) 
Rop,o = 4.13 cce./min. (9) 


Since it is known that small amounts of chloride pass into the in- 
testinal lumen under conditions closely approximating those of the 
experiment, Roo is calculated from the general fluid circuit theory. 
Using a simplified equation for the general theory (Peters, 1940, equa- 
tion 11) we have 


C(dV/dt) + V(dC/dt) =C.R; — CR,, (10) 


where C, is the effective chloride concentration of the fluid passing 
into the intestinal lumen at rate R;. dV/dt, C., R;, and R, are here 
assumed to be constants. 

If dV/dt = —D, then, from equation (10) 


H, C. PETERS 151 


@C MOR, — CR, ++ DC 


dt V Cee 
Substituting dV/—D for dt and rearranging, 


dC dV 


Cle Se IIE (12) 


Since the solution placed in the intestine was isotonic and the experi- 

ment lasted only ten minutes, osmosis is assumed to be negligible and 

DShe- hi, = kR,— D. (18) 

Substituting R, — D for R; in equation (12) we obtain 
dC E AV hed V 


Ca 3 ee eee 
Integrating between the limits C, and C and V, and V , 
Cia CoV Ve 
C= tyV,=(r,) te 


Roci can be calculated from this equation only if a value is as- 
signed to C.. The author (1941) has recently estimated upper limits 
for C. from concentration-time curves for the accumulation of chlo- 
ride in originally chloride free solutions placed in the lower ileum. 
An average upper limit of 0.018 per cent NaCl and a maximum upper 
limit of 0.033 per cent were obtained in a series of eight experiments. 
Assuming that C. = 0.018 per cent NaCl and applying equation (15) 
to the data of the chloride-bromide-deuterium oxide experiment we 
find that 


ee Coy Mins (16) 

From equations (8), (9), and (16) we obtain the relations 
Rop.0 = 1.92 Roc, (17) 
lige, ces Neiy liga (18) 


Ingraham (1935) has shown that bromide and chloride are ac- 
tively absorbed at the same rate, as would be expected from the fluid 
circuit theory. Since chloride is moving against a gradient in the pres- 
ent experiment, this may be expressed 


gr = Roc . (19) 


From equation (6) we have 
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der = Dp er Uieyon (20) 
Substituting for Roz, and dg, from equations (18) and (19), 
Opa OA coat: (21) 


In three experiments with rather porous collodion membranes at 
37.5°, dp,o/de, ratios of 1.6, 1.7, and 1.9 were obtained. Assuming 
that the diffusion areas of the intestinal epithelium have the same 
permeability ratios as the average collodion membrane used, we have 


een aes (22) 

and from equation (21) 
dp,o = 1.7 X .54Roc1 = .92Rocr. (23) 

From equation (7) 
Qp,0 = Rov,o — p,0- (24) 


Substituting from equations (17) and (23), 
p.0 — (1.92 nor 92) Root — 1.00 Rocr 0 


This is exactly what would be expected from the fluid circuit 
theory. If we assign to C, the value for the maximum upper limit, 
0.033 per cent NaCl, we obtain in the same way @p,o = 1.04R oc, and 
if C. =0, we find @p,0 = 1.00Ro:. If Ce. is assumed to have the aver- 
age value for the upper limit of 0.018 per cent and dp,o/ds, is varied 
from 1.52 to 1.88, @p,o does not deviate from the theoretical value by 
more than-ten per cent. 
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On the basis of thermodynamic considerations, expressions for 
the relation between the 7-th and j-th ions and solvent separated into 
subsystems by a semi-permeable membrane are derived. From this 
basic equation expressions for the transfer of solvent attending so- 
lute changes are developed. The application of these expressions to 
biological systems involving solvent transfer, as in cellular growth, 
chloride shift in the red blood cell, imbibition of water into cells 
is discussed. 


I. Introduction. In the usual thermodynamic derivation of the 
membrane equilibrium equation, it is tacitly assumed that no solvent 
is transferred across the membrane. (See for instance, Bolam, 1932.) 1 
The fallacy of this assumption is apparent. In fact, the necessary 
transfer of water for maintaining osmotic equilibrium is explicitly 
stated in one of the classic applications of this theory to water and 
electrolyte balance between cells and plasma in blood (Van Slyke, 
1923, and Henderson, 1928). Experimental determinations of the 
ratios, [Cl-]-/[Cl];, and [HCO;]./[HCO;], have failed to satisfy 
the equality demanded by the simple equation (Van Slyke, 1923, and 
Henderson, 1928). In the present note there is first derived a basic 
relation between the 2-th and j-th ions and the solvent, and then its 
implications of solvent transfer attending ion changes are studied 
(Donnan and Guggenheim, 1932) .? 


II. Derivation. We will begin by writing the chemical potential 
(partial molar free energy) of the i-th substance as a series of terms 
in which (at constant temperature) each variable term represents the 
energy above the standard energy possessed by 7 because of some 


1 This statement does not imply that the question of water transfer across 
membranes has been altogether neglected. See for example the important papers 
of Longsworth, 1933, 1934; Schreinemakers, 19334, 1933b, and Liu, 1931 and 
1932. 

2The basic equation for the sum of positive and negative ions has already 
been given in this beautiful paper by Donnan and Guggenheim. However, their 
equation cannot be used to derive expressions for solvent transfer attending 
changes in particular electrolytes. 
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type of work which has been performed upon it, e.g., osmotic, elec- 
trical, etc. For the system we are interested in: 


i= i + RT loga;+aVit2Fy. (1) 


This symbolism requires some comment. j;* is the standard chemical 
potential—the partial molar free energy when 4 =1,2—0, and 
y =0. a= the pressure acting throughout the phase; in particular 
this factor contains the osmotic term. In the usual derivation it is 
omitted, implying the solvent cannot move from phase to phase. 2: 
is the valence of the i-th ion. F is the Faraday constant, and y is the 


electrostatic potential of the phase. V; is the partial molar volume of 


the i-th ion, Vi = (0V/0ni)r,2n,, 7 #7 where V is the volume of the 
system and n; is the number of moles of the i-th substance. Clearly 
V; will depend on what substances are mixed with 7 in each particu- 


lar concentration, and measurements of V; must be made on a solu- 
tion similar to the one being considered by the theory; its value is 
not merely molecular weight/density.? There are theoretical difficul- 


ties involved in obtaining V; when 7 is an ion. This matter will be 
considered later (p. 155). a; is the activity (referred to mole-frac- 
tions) of the 7-th ion. Strictly speaking it is not possible, by thermo- 
dynamics alone, to obtain a; for ionic substances, and the Lewis con- 
cept of mean ionic activities (Lewis and Randall, 1923) should be 
used. However, the Debye-Hiickel theory does make possible the cal- 
culation; so that we will retain this form, bearing in mind the way 
in which a; is to be obtained. 

Using the form (1), Donnan and Guggenheim, 1932, have given 
the relation between solvent and i-th ion activities and, A y the mem- 
brane potential, or electrical potential difference between the two so- 
lutions ; 


1 an RT re : 
= RT log = — log : +S ay, (2) 


, ff 


al 1 a a; V 


Manipulating expressions for the j-th ion in a similar way and elimi- 
nating A y between the two equations of the type (2), we obtain: 


£0 (Vic \ waters aera ome 
as [ee Oe = log ee ee 2 
A Cd rea a 


It is important to note that the first term of equation (3) represents a 
correction, as it were, on the simple membrane equilibrium equation; 


8It is unfortunate that this concept has been confused in otherwise good 
expositions. See, for instance, the calculation of Liu (1931) and of Adair (1937). 
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a" ;/U; = 0"5/a;. (4) 


III, Deductions. Ion Ratios: It is clear that equation (3) redutes 
to the simple form in the case that V;/z;= Viz; onit ds — 0, ein 
systems consisting of dilute solutions a’, ~ a',; however, systems 
wherein one phase contains non-permeant ions in fair amounts may 
present considerable deviations from this equality. Again, whenever 
a marked concentration of substances occurs, the ratio, @";/a’,, will 
depart from unity and consequent deviations from the equilibrium 
ratios predicted by the simple equation would be expected. Cases in 
point are the well-known Valonia cells, Ascidian organisms, red blood 
cells, ete. 

Numerical estimates of the term, 1/a where 


11 %_Y, 2 
demand a knowledge of the partial molar volumes of ions. Partial mo- 
lar volumes of individual ions exhibit the same thermodynamic inde- 
terminacy as activities of individual ions. As in studying the latter, 
recourse to molecular theory may be helpful. Partial molar volumes 
for electrolytes of biological importance, e.g., NaCl, KCl, ete. 
are of course readily available in the literature. We can proceed by 
postulating that for some electrolyte the V of the positive ion may 
bear some relation to the V of the negative ion. On the basis of the 
striking symmetry of KCl (e.g., equal migration velocity of the two 
ions in water, geometrical symmetry of KCl crystals, etc.) we as- 
sumed that in KCl, 
View = Vei- (6) 
and 
Vent Vala Known. (7) 


One can then find the V’s of other ions by the simple procedure of 
solving an equation of the type of (7). For instance, 
Lo oe Ue eee eae Vat ae Ve : 

In this manner a table of “ionic partial molar volumes” can be built 
up. It may be remarked that cursory application of this procedure 
led to reasonable agreement between values for V; calculated from the 
V of different compounds containing the 7-th ion. 

“Chloride shift’: It is found experimentally that during the oxy- 


4 We are indebted to Dr. W. F. Libby for valuable suggestions. 


156 MATHEMATICAL BIOPHYSICS 


genation and reduction of hemoglobin in red blood cells, transfer of 
bicarbonate and chloride ions across the cell membrane occurs in op- 
posite directions. This phenomenon is commonly known as the “chlo- 
ride shift’’.. Under these circumstances, water is simultaneously trans- 
ferred. Despite this evidence of osmotic work, the simple equation 
[Cl]./[Cl], = [BHCO,]-/[BHCO,],, as derived on the assumption 
that no osmotic work is done, has been traditionally used as an ex- 
planatory principle. Substitution in it shows quite obviously that 
either of the ion ratios varies directly with the other; i.e., Cl and 
HCO, should transfer in the same direction, which is contrary to fact. 
However, if the general equation (3) is used, no such antinomy is 
reached, since the transfer of solvent can account for the passage of 
the two ions (Cl and HCO;) in opposite directions. 

Solvent Transfer: A Suggested Mechanism: The uptake and ex- 
trusion of water, e.g., water uptake in growth, etc., is common to all 
metabolizing cells. The mechanisms of this process may be various, 
but certainly equation, (3) clearly indicates that changes in solute 
concentrations effect water transfers. (It should be noted that small 
changes in the activities of water involve relatively large amounts of 
water transferred, for the reason that water exists at a much greater 
concentration than any solute.) In this section we will derive simple 
expressions for the amounts of solvent transferred when the solute 
concentrations change in certain ways. All these expressions are for 
small variations in solute concentration. The method for obtaining 
the corresponding expressions for large variations is obvious. Let us 
suppose that f is the water concentration in a cell, and that it can be 
written as a function of two concentrations of solutes x and y. Thus: 


fST(e, yn (8) 


The increment in f attending small variations in x and y, is 


sf=(Sa)o0 +(sz) oz. (9) 


The present treatment is made on the assumption that the cell is sus- 
pended in media which maintain constant concentration (infinite res- 
ervoirs). Double primes denote cell interior; single primes, the me- 
dium in which the cell exists. Furthermore, for the sake of formal 
simplicity, we will now assume ideal solutions, whence we may re- 
place a; by N;. 

A. Solvent Transfer with Ion Changes. (i) Solvent and Ion 


transfers: Rearranging equation (3) preparatory to differentiation 
and substituting from equation (5), we have 
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log Nx" — log Ni’ = (a/2:) log (Ni"/Ni’) — (a/zi) log (N;"/N;') . 
The increment in N,” is then given by 

6 (log Nx")=(6 Nx"/6 Ni")=(a/zi) (6 Ni"/Ni") — (a/2;) (6N;"/N"5). 
(10) 
Equation (10) is quite general, and gives the per cent changes in sol- 
vent mole fractions when the mole fractions of the i-th and j-th ions 
undergo variations, 6 N;” and 6 N,;” respectively. These variations 
are not, however, completely independent. For electrical neutrality, 
40n, + 2;6n;=—0, CLE) 


where of course ; and n; represent the numbers of moles of the i-th 
and j-th substances. Now since the total number of moles, A , remains 
fairly constant we may divide equation 11 by A and obtain 


z0N,+2;6N;20. (12) 
From equation (12) it follows that 6 N; = — - 6 Ni. 
j 
Substituting into equation (10) we obtain 
6 N,’ 1 Zi 
SS —s ———_ N;". 1 
N,” 8 INE ¥ Cie wr)? ( 2) 


Rearranging within the brackets of equation (13), and substituting 
within it the explicit value of a, 


” V Bi Nua ae IN, 
6N, nite Vi F J ) Nee (14) 


NGO zj N,;" N,;” Zj V; = Zi V; 


If we change the i-th and the j-th mole fractions and the mole frac- 
tions of no other permeant ions or non-ionized substances, the per 
cent change in the mole fraction of solvent is given by the above equa- 
tion (14). It might be well to pause here and consider how such a 
change may be brought about. Returning to equation (11) we see 
that z; and z; are of the same sign, n; and n; must be of opposite signs ; 
in other words, the exchange of ions with the external medium can 
lead to solvent transfer. On the other hand, if z; and z; are of oppo- 
site signs, then n; and n; must be of the same sign, implying that the 
transfer in toto of an electrolyte also effects a solvent transfer. It 
might at first sight be thought that the apparently arbitrary changes 
we have just mentioned violate the mole fraction equation, 


Sy NZ=1 (15) 
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This is not the case, however, provided we realize that in equation 
(15) the summation must be made over all substances, and not only 
over the permeant ones. Into the expressions of the Donnan Equilib- 
rium only the permeant substances enter. We might re-write equa- 
tion (15) as follows: 
> N,(impermeant) + 5 Na(permeant) + Ni; +N;=1. (16) 
uw d#i 


dzj 


In equation (16), N; and N; are changed according to our discus- 
sions; Sz must remain constant by hypothesis, but }, can compensate 
the change in the i-th and the j-th substances and thus preserve the 
equality. 

(ii) Changes in the Membrane Potential Attending Solvent and 
Ion Transfers: Since by rearranging equation (2) we have, 
Vik T/L ore NG eee, 

( lo ) 


es — — log — 
2,F Ke Nx V; N;’ 


it follows that 


VR (ae Nal tgs laa 
s(Ay) = (= 55 ) 17 
" slailica ATEN ViaieN teens Vie oe 
Substituting from equation (14) into equation (17), 
| ViRT (1 (2P@Nj" +22N," 
d(Ay) =o 4 =) (18) 
aN," F 2; N; 2; Vi — & V; Va 


Equation (18) now gives the change in membrane potential when 
N; and N; vary subject to the conditions discussed above. It is clear 
that in general the changes occurring in our systems can be described 
by any two of the three functions: solvent transfer, ion transfer, mem- 
brane potential. 

B. Solvent Transfer Attending Concentration Changes of Undis- 
sociated Substances. Donnan and Guggenheim, 1932, give the rela- 
tion between the concentration of the w-th undissociated substance 
and the solvent as, 


1 iri 1 1 N;° 
—a —— 116, rh 
Va N Go ave N, oe 
Taking increments as before, 
O. Nac he Vets Ne" 
(20) 


Pe ————_, 
N,’ Ve Ne 
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We can imagine that the mole fraction of the u-th substance may 
change either by uptake or by extrusion of the u-th substance, but 
perhaps biologically more important is the possible polymerization of 
the u-th substance into an impermeant material according to the re- 
action 


nm Ue ereenet) ae ner cant) ° (21) 
Implying of course, that 
ON," =moéN,,. (22) 


If m is very large, the effect of this mechanism of solvent trans- 
fer would be correspondingly marked. An interesting reaction which 
suggests itself is the so-called “dehydrolytic synthesis” of peptides, 
proteins, etc. In general if a peptide linked molecule contains m mono- 
equivalent amino acids, (m — 1) moles of water were produced by its 
polymeric synthesis. Thus, if 6 N,” is the increment in impermeant 
ion, the maximum increment in amino acid would be m 6 N,”. (In 
the case where only one acid is involved), and the increment in water 
is (m—1) 6 N,” =6 N;,” (reaction). Now in equation (20) , 


%) N,” =0 Ne Geeeetien) + 6 Ny" (transter) = 
So we may write: 


ON," (transfer) _ V8 N,” PPROUN oe 


- (reaction) (23) 
N,” Vu N,," N, 
multiplying and dividing the last term in equation (23) by 6 N,", 
and remembering N,” (reaction) /N,,” = m—1/m, we have 
6 Ny" (transfer) _ : Vitek tee =) 5N,". (24) 
Nae Vie oe et NGL 


In cases where the molecule is being polymerized from more than one 
acid, similar though more awkward expressions than equation (24) 
could be at once written down merely by breaking up the N,” into 
fractions. 
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BOOK REVIEWS 


HANS SCHAEFER. Elektrophysiologie. I. Band. Allgemeine Elektrophysiologie. Vi- 
enna: Franz Deuticke, 1940, XI +. 522 pp. 94 figs. 

BERNHARD KATZ. Electric Excitation of Nerve. London: Oxford University Press, 
1939, IX + 151 pp. 31 figs. 


An unsuspecting reader opening either of these excellent volumes at random 
would very likely think that he had before him a physics monograph rather than 
a physiology text. And, indeed, this is a very encouraging state of affairs, for 
each of these books abundantly proves that electrophysiology has finally and ir- 
revocably adopted the method, if not the aim, of physics, both experimental and 
theoretical—so much so that electrophysiology is certainly one borderline field in 
which philosophic distinctions between biology and physics have been ignored. 

The scope of Dr. Katz’s book is considerably smaller than that of Prof. 
Schaefer, and, as a result, his exposition is somewhat more unified than the lat- 
ter’s. Beginning with some fundamental physiological observations (local summa- 
tion, ete.), Dr. Katz proceeds to a formulation of the current two-factor excita- 
tion equations together with a discussion of their experimental support. The 
viewpoint is formal throughout, no particular assumption being made concerning 
the physical nature of the excitation parameters. In the second half of his book, 
Dr. Katz considers those aspects of nerve activity, such as electrotonus, conduc- 
tion, local response, ete., which are manifestations of the kernleiter structure of 
the nerve. The viewpoint here can no longer be strictly formal, and Dr. Katz’s 
attempt to incorporate the formal, linear excitation equations into a description of 
sub-liminal response leads to certain logical inconsistencies. Either one must 
follow Rushton and work out a completely physical theory which includes a uni- 
fied discussion of both excitation and conduction, or one must ignore the kernleiter 
properties of the nerve and treat subliminal as well as superliminal responses on 
a purely formal phenomenological basis. The latter course leads to non-linear 
excitation equations, which, though complicated, seem to afford the possibility of 
including the phenomena of repetition and recovery into the mathematical de- 
scription of nerve activity. 

Prof. Schaefer’s treatise (it contains no less than 6000 separate references, 
all neatly classified and indexed) takes its place as the obvious successor to Bern- 
stein’s Elektrobiologie and Biedermann’s Elektrophysiologie. The work begins 
with a very complete survey of the passive-electric properties of living tissue and 
then goes on to an exposition, both mathematical and experimental, of the two- 
factor theory, and of local response. As in Dr. Katz’s book, no consideration is 
given to the necessity for a non-linear modification of the theory to account for 
events occurring near excitation. Professor Schaefer identifies the parameter of 
excitation with the membrane potential difference, and this leads him tentatively 
to associate the process of accommodation with the spread of potential away 
from the site of excitation due to the geometry of the excitable structure. Such 
an assumption concerning accommodation is certainly conjectural; its mathe- 
matical consequences (which could be tested experimentally) have not been 
worked out yet. 
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A well-organized, clear discussion of the characteristics and possible origin 
of various biological potentials (demarcation potential, action potential, mem- 
brane potential) occupies almost a third of the text. It is a little disappointing 
that the theoretical discussions of conduction which have appeared at various 
times in the literature are barely mentioned; to this reviewer it seems that they 
are comparable in importance to the various excitation theories, since tentative 
conduction equations relating velocity, safety factor, and physical constants of 
nerve could serve as a framework for experimentation on conduction of the im- 
pulse in the same way that the two-factor theory has served for organizing data 
on excitation. 

The volume ends with a consideration of synaptic transmission, and the effect 
of unphysiologic conditions on bioelectric phenomena. Professor Schaefer in- 
clines very heavily toward an electrical interpretation of synaptic transmission, 
and he presents his evidence well. However, it should be pointed out that in a 
science as young as electrophysiology (one can argue that the modern science 
started with the advent of the cathode ray oscillograph) an author is apt, un- 
consciously to be sure, to allow his own experiments and interpretations to affect 
his judgment of the relative importance of the data of others; and, to some ex- 
tent, Professor Schaefer seems to have fallen into this error whenever he con- 
tends that the excitation process “is” the recorded action potential. 

But these are very minor criticisms for a work of this magnitude. This vol- 
ume is, in the reviewer’s opinion, the most important modern book on bioelectric 
phenomena to have appeared. It is sincerely hoped that some anonymous bene- 
factor of English-speaking electrophysiologists will gain the everlasting grati- 
tude of his colleagues by translating the book into English. 

ALVIN M. WEINBERG 
The University of Chicago 
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